Atomistic study of electronic structure of PbSe nanowires 
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Lead Selenide (PbSe) is an attractive 'IV- VI' semiconductor material to design optical sensors, 
lasers and thermoelectric devices. Improved fabrication of PbSe nanowires (NWs) enables the 
utilization of low dimensional quantum effects. The effect of cross-section size (W) and channel 
orientation on the bandstructure of PbSe NWs is studied using an 18 band sp'^d^ tight-binding 
theory. The bandgap increases almost with the inverse of the W for all the orientations indicating 
a weak symmetry dependence. [Ill] and [110] NWs show higher ballistic conductance for the 
conduction and valence band compared to [100] NWs due to the significant splitting of the projected 
L-valleys in [100] NWs. 
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FIG. 1. Projected PbSe unitcell with wire axis orientations (X) along (a) [100], (b) [110] and (c) [111]. The cross-section size 
is of 2.5nm (Y) x 2.5nm (Z). Lead (Pb) and selenium (Se) atoms are shown in gold and grey color respectively. 



Appealing PbSe bulk properties: Lead selenide (PbSe) is a narrow, direct band gap semi-conductor material 
j^gbuik _ 0.16eV at 4K [D [2]) with useful electrical, optical and lattice properties [3H5]- It is used extensively in 
optical devices [31 [S], lasers [3 [5] and thermoelectric devices [H 13 En]- The large Bohr exciton radius of about 
46nm in PbSe makes it a suitable system to study quantum confinement effects on electrons and holes [6l [TTUI3] . 
Recent progress in multiple exciton generation (MEG) in PbSe with higher optical efficiency has renewed interest in 
the optical properties of PbSe [5 . Extremely low thermal conductivity of PbSe (~ 2 W/m-K in bulk PI3] to ~ 0.8 
W/m-K in NWs [TS]) also makes it a suitable thermoelectric material [HII^. PbSe can become a preferable material 
over Lead Telluride (PbTe) due to the higher availability of selenium (Se) compared to tellurium (Te) ;16]. 

Nana structures PbSe : One dimensional nanostructures of PbSe like nanowires (NW) and nanorods (NR) combine 
the interesting bulk material properties as well as the quantum confinement effects which can lead to better thermo- 
electric [ini HH US] and optical [3] devices. The analysis of the physical properties in PbSe NWs will require proper 
understanding of the electronic structure, which is the theme of the present paper. 

PbSe NW growth technologies: With recent advances in the growth and process technology, the fabrication of PbSe 
NWs have become very efficient and controlled. PbSe NWs are developed using a variety of methods like, chemical 
vapor transport (CVT) method oriented nano-particle attachment [TT], electro-deposition without catalyst [T^ 
and with catalyst [TT', hyper- branching [18 , growing PbSe structures on phosphate glass [19], etc. These methods 
enable PbSe NWs fabrication with a variety of growth directions and surfaces. 

PbSe crystal structure: Bulk PbSe has a stable rocksalt structure with a co-ordination number of six at room 
temperature (300K) and normal atmospheric pressure [I1[I1]. The lattice constant is 0.6121 nm [1' at T = 4K, which 
is also utilized in our bandstructure calculations. The PbSe NWs are constructed using the same bulk structure with 
three different wire axis orientations of [100], [110] and [111] (Fig. [ij. 

Electronic Structure calculation: The main methods for analyzing the lead salt structures are (1) continuum method 
like 4/8-band k-p ^6], (2) semi-empirical atomistic methods like tight-binding (TB)[5J[2n], pseudo-potential approaches 
[2T] or (3) first principles methods [TS] [TBI [22l [23] . The 4/8-band k-p method captures some quantum confinement 
effects, however, it fails to capture the interaction of various valleys present along the L-K and L-X directions, the 
band anisotropy and the correct frequency dependent dielectric function [131 1201 121] . The first principle calculations 
are highly accurate but the demand for computational power is very high and these methods are limited to solving 
only small structures with few thousand atoms [201 IM] ■ The semi-empirical methods are versatile in terms of the 
involved physics and can handle a large number of atoms (> 10 million atoms [24] ') making them suitable for electronic 
structure and transport calculation in realistic device structures. However, an integral part of these methods is the 
requirement of correct semi-empirical parameter sets to properly represent the electronic structure properties like 
bandgaps, effective masses, wave-function symmetry, etc. In this work we utilize a semi-empirical atomistic TB 
method based on sp^(f formulation with spin orbit coupling (SOC) optimized for bulk Pb salts by Lent et al [5] to 
calculate the electronic structure in NWs. The inclusion of SOC is important in PbSe since both CB and VB have 
strong p-contributions from Pb and Se atoms, respectively [T31I20] . 

Surface atoms: The TB parameters developed for bulk Pb salts [2] are also applied to nanostructures with finite 
boundaries. The surface atoms are not passivated for Pb salt nanostructures. It has been previously shown that 
surfaces do not introduce states in the energy gaps of the bulk band structure, even though they are not passivated 
[20] and the same is obtained in the present study. The surface states are mainly p-like for Pb salts which are strongly 
coupled with orbitals of atoms at the interior unlike the zinc-blende semiconductors where atoms are described by 
hybrid sp^ orbitals which remain uncoupled at surfaces, forming dangling bonds [25] . The lack of surface passivation 
has also been pointed out by the first principle calculations in stoichiometric Pb salt nanostrutures [131 123| . 
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FIG. 2. Bandstructure of PbSe NWs for (a) [100], (b) [110] and (c) [111] oriented channels. The CB minima and VB maxima 
are normalized to zero for simplified band comparison. Only the first 5-10 sub-bands are shown. 



Electronic Structure: The bulls; effective masses calculated using the TB parametrization in Ref. 2] is obtained for 
(i) electrons as mU/m^ = 0.087(0.07)/0.036(0.04) = 2.4 (1.85) and (ii) holes as m|[/m^ = 0.094(0.068)/0.031(0.34) = 
2.9 (2.0). The values in the parenthesis are from Ref.Jl]. The TB parametrization in Ref. ^{TBa) is preferred over 
the parametrization in Ref. |20) (TBb) since the TBj^ model captures the mass anisotropy quite well, an important 
requirement for TB models as pointed in Ref. pT] . 

The bandstructure of conduction (CB) and valence band (VB) for PbSe nanowires are shown in Fig. [2] for 3 
different wire orientations. For all the wires the conduction band minima (CBM) and valence band maxima (VBM) 
are normalized to zero to enable a better comparison of the valleys. Few important points to observe are, (i) bulk 
L valleys are projected at the Brillioun zone (BZ) edge at X for all orientations (Fig. [2]). The [110] wire has an 
additional projected valley at the F position (Figj2]b, e), (ii) [100] NW show a degeneracy (D) of 1 for both CB and 
VB due to strong valley splitting, however [110] and [111] NWs show a D of nearly 2 (extremely small valley splitting) 
for the positive 'k' states. 

Bandgap and bandedge variation: The CBM and VBM variation with cross-section size (W) and orientation are 
shown in Fig. |3] a. All the orientations show quite similar band-edge variation with W. As the cross-section size 
decreases the geometrical confinement increases which pushes the CBM (VBM) higher (lower) in energy. The variation 
in the bandedges with W also compare surprisingly well with a previous 4-band k-p calculation done for [111] cylindrical 
PbSe NWs [H] . The bandgap variation with W in PbSe NWs can be fitted to the following analytical expression. 



Eg^P/iWr, 
where a represents the power law dependence on W. 
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FIG. 3. (a) Bandedges for square PbSe nanowires with [100], [110], and [111] channel orientations. Bandedge result using 4- 
band k-p calculation for cylindrical PbSe nanowires from Ref. [26]. (b) Bandgap variation for all the nanowires. Dots represent 
calculated TB values and lines represent analytical fitting with cross-section size. 
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FIG. 4. Ballistic conductance in PbSe NWs for (a) CB and (b) VB for 3 different wire orientations at T = lOK. The temperature 
comes from the Fermi-Dirac distribution of the carriers. 



For [111] NWs the value of (3, {a) is 1.726 (0.8734). For [110] NWs these values are (a)) 1.564 (0.8592) and for 
[100] NWs these values are (/3, (a)) 1.87(1.072). The bandgap values (Fig. |3]b) roughly show an inverse relation with 
W for all the NW orientations which is very different from the prediction of effective mass approximation (Eg oc W~'^). 
Similar results for Eg have been obtained by other independent calculations carried out in PbSe nanostructures using 
first principle calculations [13] as well as TB calculations [2D]. This justifies the application of TB electronic structure 
calculation which correctly captures the quantum confinement effects in ultra-scaled PbSe NWs. 

Ballistic conductance: Transport properties of PbSe NWs are revealed by the electronic conductance which is 
calculated using Landauer's formula [27|. Figure H shows the normalized ID ballistic conductance for electrons in the 
PbSe NWs for 3 different orientations. [110] and[lll] oriented wires show higher G value for both the CB and the 
VB compared to [100] NWs. A larger valley splitting in the CB and the VB in [100] NWs decrease the conductance 
compared to the other two orientations. The normalized conductance increases in steps of 2 for [100] NWs for the 
CB but not for the VB (Fig. |4]) which shows that the CB and the VB are not entirely symmetric in energy, a result 
similar to the one given in Ref. |21j . Thus, transport characteristics show the influence of geometrical confinement 
and channel orientation. 

Summary and Outlook: PbSe NWs have tremendous potential of becoming the next generation thermoelectric 
and optical devices. The proper understanding of the physical properties of these ultra-scaled PbSe nanowires will 
depend strongly on the correct electronic structure calculation. We have presented the application of semi-empirical 
tight-binding theory to these NWs to understand the variation of conduction and valence bands and the position of 
the important energy valleys. The variation of the bandgap and the bandedges with NW cross-section size is strongly 
influenced by the cross-section size and channel orientation. Simple effective mass theory cannot predict the bandgap 
variation. The amount of valley splitting is strongly dependent on the type of geometrical confinement as reflected 
in the ballistic conductance of the PbSe NWs. [110] NWs provide the maximum ballistic conductance for both the 
CB and VB. The tight-binding analysis of electronic structure opens door to explore the optical and thermoelectric 
effects in PbSe NWs. 
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